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ABSTRACT 



Boiling Heat Transfer with Low Steam Quality 
Thomas Frank Karucci 

Submitted to the Department of Naval Architecture and Marine Engineer- 
ing on Kay 21, 1965 , in partial fullfillment of the requirements for 
the Master of Science Degree in Mechanical Engineering and the Profes- 
sional Degree, Naval Engineer. 



Heat transfer for water and steam mixtures in a horizontal stain- 
less steel tube under non-boiling, surface-boiling, and fully-developed 
boiling conditions was investigated. 

Local heat transfer coefficients are presented for a limited range 
of flow parameters. Inside diameter of 0.21*25 inch and a range of Ih.h 
to 78.5 L/D are considered. Exit pressure varied from 60 to about 72 
psia for mass velocities from 5.9 x 10^, to 2.38 x 10^ lbs/hr ft^. The 
range of steam quality covered was 0 to hh%» The majority of this data 
was taken at inlet temperatures of about 290°F. Heat fluxes ranged from 
2 x 10 ? to 1.8 x I06 Btu/hr ft2 and were limited by burnout conditions* 

Heat transfer results for non-boiling and subcooled nucleate boil- 
ing are in agreement with previous investigations. Bulk boiling data 
is presented at the two lower flow rates. This data shows that as quali- 
ty is increased, the local heat transfer coefficient also increases. This 
is indicated on the boiling curve as a shift toi-rard loxcer wall superheat 
as quality increases. Similar results are indicated to a lesser degree 
with subcooled liquids at low values of subcooling. This appears to be 
due to the relatively large volume fraction of vapor created in the nu- 
cleate boiling process. 

Experimental heat transfer coefficients obtained are considered to 
be in fair agreement with predicted values. 
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Title: Assistant Professor of Mechanical Engineering 
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NOMENCLATURE 


A 


c surface area 


A c 


■ test section cross section area 


C P 


1 specific heat 


D 


1 diameter 


E 


1 voltage across the test section 


E s 


1 shunt voltage 


G 


1 mass velocity 


h* 


enthalpy 


h 


heat transfer coefficient 


I 


current 


k 


thermal conductivity 


L 


heated length 


P 


pressure 


q 


heat transfer rate 


q/A 


heat flux 


(q/A)i 


heat flux at the inception of boiling 


(q/A) cr - 


critical heat flux 


(qA) F c 


forced-convection heat flux 


R 


electrical resistance 
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T 


IK 


temperature 


T 

1 w 


BK 


inside wall temperature 


T 

ov 


- 


outside wall temperature 


V 


a 


velocity 


w 


- 


mass rate of flow 


X' 


ss 


position along the test section measured 
inlet to the heated section 


X 


m 


steam quality 


AT 


m 


temperature drop through the tube wall 


A "^sat 


S3 


wall - minus - saturation temperature 


A -^sub 


m 


saturation - minus - bulk temperature 


M 


m 


dynamic viscosity 


cr 


IK 


surface tension 


P 


K 


density 


Pe 


SC 


resistivity 


X 


m 


latent heat of vaporization 






Dimensionless Groups 


Nu 


a 


Nusselt number *■ hD/k 


Pr 


SB 


Prandtl number * C 0 M/k 


Re 


m 


Reynolds number “ &T>/M 



Subscripts 



b 

BO 

i 

exit 

L 

sat 

t 

v 

w 



bulk or mixed-mean conditions 

burnout condition 

inlet bulk condition 

test section exit bulk condition 

liquid condition 

saturation condition 

tube 

vapor condition 

condition at heat transfer surface 



INTRODUCTION 



The current development of high performance heat-transfer 
equipment has stimulated a great amount of research in boiling 
heat transfer. Past studies of boiling heat transfer have been 
a factor in the development of nuclear reactors and rocket motors. 
Vaporization of the coolant is frequently desired, such as in the 
boiling-water reactor and bootstrap rocket engine. 

There has been a vast amount of work in the field on forced 
convection-boiling heat transfer. Kost of this has been concen- 
trated on the prediction of burnout under given conditions of 
fluid, surface, and flow parameters. These studies have produced 
numerous correlations; the best of which can predict critical heat 
flux to about plus or minus 30$ of experimental data. 

Within the past ten years there has been an increase in the 
number of investigations designed to study heat transfer rates en- 
countered in boiling heat transfer. The results of these studies 
have been presented in two basic forms, ^irst, in the form of the 
boiling curve where heat flux is plotted versus wall superheat 
(T w -T sa t)» Th; i s has usually been utilized where the bulk fluid is 
subcooled, and indicates that heat flux can be represented as: 

^ - C (T„ - T sat ) n (1) 



Other studies considered bulk boiling and present results that show 
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the local heat transfer coefficient (h) increasing as steam quality 
increases beyond a certain range. This range has been generally 
agreed upon as beginning at about l$% quality. As quality varies 
below this value h is constant. 

Generally, when liquid flows through a heated tube, subcooled 
nucleate boiling can take place before the bulk of the fluid reaches 
its saturation point, if the impressed heat flux is high enough. The 
equilibrium thermodynamic quality would be zero in this subcooled 
region of flow, but the void fraction could be substantial, depending 
of course on the fluid properties and heat transfer rate. As flow 
proceeds down the tube the bulk of the fluid becomes saturated, net 
vapor generation commences. In this portion of the flow, the vapor 
generation can take place by additional nucleate boiling and/or by 
evaporation at the liquid-vapor interface. If the tube is long enough, 
it is possible to reach the point of 100$ vapor. 

This study was conducted to extend the fully-developed boiling 
curve in the region of bulk boiling and to determine what effect of 
variation in h with quality has on the boiling curve. It was also 
pointed out by Professor 3ergles that non equilibrium voids created 
in subcooled region should have an effect on h. 

The equipment used limits the study to low pressure, low flow 
rates and moderate heat fluxes. Fluid conditions included low sub- 
cooling and low steam quality. This study is particularly applicable 
in the development of evaporators for saline water conversion. 
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PROCKDURE 

Operation and Data Taking 

The experimental apparatus (Figure 1) used in this study was 
one already in operation at the M.I.T. Heat Transfer Laboratory. 
Only the test sections had to be fabricated to suit this study. 
Generally, the test sections had to be instrumented to provide 
means of determining local pressure and the temperature at the 
heat transfer surface. To obtain these measurements, pressure 
taps were located along the tube and thermocouples were placed 
on the outside wall of the test section. To insure that the out- 
side wall was adiabatic, two guard heaters were used and main- 
tained at a temperature equal to the tube wall temperature. A 
brief discussion of the experimental apparatus and a detailed 
discription of the construction of the test section is presented 
in Appendix A. A diagram of the heat transfer test section is 
shown in Figure 2. A preliminary investigation of critical heat 
flux was conducted to be able to predict the critical heat flux 
for the range of parameters to be studied. This was done pri- 
marily to avoid burnout of the elaborately instrumented heat 
transfer test sections. The test sections used to determine 
critical heat flux were not instrumented to obtain wall temp- 
eratures and were constructed with only one pressure tap located 
at the test section exit. 

Prior to installation, the test section was given a final 
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internal cleaning with acetone. After installation, the entire test 
section was wrapped with insulation to minimize external heat losses. 

The loop was then filled with distilled water and all lines vented to 
remove air. The degassing operation was then initiated by heating the 
water in the degassing tank. With the water boiling, the loop water 
was heated and slowly sprayed into the top of the tank. Results from 
Reference (1) indicated that 20 minutes of such degassing reduced the 
dissolved-air content to approximately 1 c.c. /liter. Two phase flow 
can influence system stability and recent studies have indicated that 
system instability can cause premature burnout. Air trapped in the 
system would provide a certain amount of compressibility. Reference (2) 
discusses system instabilities in more detail. After degassing, the 
heat exchanger cooling water was turned on to bring the system to nor- 
mal temperature. The loop was then ready for operation, circulating 
distilled, deionized, and degassed water. 

Critical heat flux daoa was obtained by increasing power gradu- 
ally while maintaining constant flow rate, exit pressure, and inlet 
temperature. Critical heat flux was evaluated at the last recorded 
conditions prior to physical burnout of the test section. To vary 
exit quality conditions, inlet temperature was varied in subsequent 
runs. 

The general procedure for taking heat transfer data was similar 
to that used for critical heat flux data. The major difference was that 
at each power level, while other conditions were maintained, the guard- 
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heater power was adjusted to the balance point, i.e., where the average 
guard-heater temperature (two- thermocouples) was the same as the aver- 
age outer tube wall temperature (three thermocouples). These tempera- 
tures could not be kept exactly equal, but were maintained to within 
plus or minus 5°F of each other. At equilibrium, wall temperatures, 
fluid temperatures, pressures, pressure differences and test section 
voltage and current were recorded. The heat flux was then increased 
to the next level and the process repeated. 

After the highest heat flux was completed, the power was reduced 
and the inlet temperature changed. Readings were then taken for a set 
of heat fluxes at a new inlet temperature. 

The maximum attainable heat flux for any given run was determined 
by burnout except in one instance due to voltage limitations of the 
generators. The burnout heat flux was estimated from data presented 
in Reference (3) and from data taken in this study. It was desired 
to approach but not reach burnout, temperature and pressure levels 
were largely determined by the system. Beginning at the lowest inlet 
temperature attainable, the inlet temperature was increased to about 
20°F below the saturation temperature for the corresponding pressure. 

The exit pressure was maintained at 60 psia except at final runs where 
the minimum exit pressure obtainable for the desired flow rate exceeded 
60 psia. This was attributed to the fact that choked flow conditions 
were present due to the large volume of vapor content in the test section. 

Initial data taken with T.S.-l and 2 was at relatively low values 
of heat flux and at low inlet temperatures. This was to obtain non- 






boiling heat transfer data which is well correlated and hence serves 
as a check on instrumentation. By increasing the heat flux slowly it 
was possible to obtain data on the inception of nucleate boiling and to 
be able to trace the boiling curve in the transition region between non- 
boiling and fully-developed subcooled nucleate boiling. 

At the higher heat fluxes with intermediate inlet temperatures the 
fully-developed subcooled nucleate boiling curve could be well estab- 
lished. 

The final data was then taken with T.S.-2 at inlet temperatures 
near saturation. The inlet temperature was maintained below saturation 
at all times since the equipment did not allow determination of inlet 
conditions when net vapor existed at the inlet. This procedure pro- 
vided data for bulk boiling over the quality range from 0 to 2k%. To 
extend the data to higher steam qualities it was necessary to fabri- 
cate T.S.-3. T.S.-3 was constructed similar to T.S.-l and 2 except 

that it was twice as long. The location of pressure taps and wall 
thermocouples for each test section is given in Appendix A. Due to 
inconsistent results as discussed later, T.S.-3 was modified by moving 
the thermocouples at the upstream locations to positions corresponding 
to L/D's at the downstream position of T.S.-2. One series of runs was 
conducted with this test section at G = 1.19 x 10^ and T^ ^ 290°F o 
Accidental burnout and time limitations precluded going further with 
this study. 
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Data Reduction 

All data reduction was done by hand calculation with the aid of a 
desk calculator. This included computation of local flow and heat trans- 
fer parameters, including flow rate, heat flux, wall superheat, bulk temp- 
erature or enthalpy as well as the dimensionless groupings Re, Nu, Pr, 
Nu/?r 0 °^. The mass velocity (G) was determined from the fluid flow rate 
and tube dimensions. The temperatures were converted from thermocouple 
emfs with NBS standard calibration. The heat input was obtained from 
the tube voltage and current (SI) and checked by the resistance and cur- 
rent, a third value was obtained from the mass flow, specific heat, and 
temperature rise. Values obtained by each were in excellent agreement 
and differed at most by 10$, but were usually within 5$ of each other. 

The latter method was of no value when the exit temperature reached sat- 
uration. The average heat flux was calculated using the several values 
of heat input. 

Heat transfer results were based on measurements of the outer wall 
temperature. A correction for temperature drop through the tube was made 
to first obtain the inner wall temperature. A modified form of the Krieth 
and Summerfield solution was used to obtain the necessary correction. 

This solution assumes an infinite cylindrical resistor, with heat trans- 
fer at the inner surface and an adiabatic outer surface. The final form 
of the equations as well as thermal-property information is given in 
Reference (1)* The final form of the curve (^T w versus I at various values 
of T^ is shown in (Figure 3). It should be noted that at high values 
of current, the wall temperature drop is extremely large. The data taken 
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in the last run (Table IV Series N Appendix C.) points out the effects 
of this on the results. The wall temperature drop is so large that the 
inside wall temperature is actually reduced as the power level is in- 
creased. Results then indicate a substantial reduction in the degree 
of wall superheat. This leads to the conclusion that the equation for 
wall temperature drop may not produce accurate results at high values 
of current. High currents were reached with T.S.-2 as well, (Table IV 
Series F) but with a shorter test section the outside wall temperatures 
were also higher, hence, the effect was not as pronounced. At each 
thermocouple position, calculations were made to obtain the local heat- 
transfer results. The fluid bulk temperatures are based on a linear 
temperature gradient where data is for subcooled bulk fluid. Where net 
vapor existed, the bulk temperature is the saturation temperature for the 
particular local pressure. The local pressure was obtained graphically 
by drawing a smooth curve through the observed pressures at three loca- 
tions. All fluid properties are based on these local bulk temperatures. 
With the saturation temperature (^ sa ^), the degree of subcooling (AT g 
T s at ~ T^) and the wall superheat ( AT sa ^ ** T^ - T g _ J .) were available. 

A sample set of calculations is shown in Appendix B. A tabulated 
summary of the data obtained is given in Appendix C„ 
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RESULTS 

Critical Heat Flux 

Burnout data were required to set an upper limit on the heat flux 

that could safely be applied to the heat transfer section. This data 

is presented graphically in Figure li, for G - 1.19 x 10^ lb/hr ft 2 

(V c^Sf t/sec). Also plotted are data taken by Wessel (3)* for the sane 

tube size with varying fluid velocity and a slightly larger tube with 

varying L/D(V ^ 5it/sec). The data shows an increase in (q/A) as 

o r 

exit quality increases. 

Heat Transfer Results 

The non-boiling results of this study are presented in Figure 5 > 
along with the line indicating McAdams (lj) correlating equation. Data 
in Figure U were obtained for G = 2.38, 1.19 and 0.59 x 10^ lb/hr ft 2 
(Vc*10, 5 > and 2.5 ft/sec). The span of Reynolds number also results 
from the bulk temperature variation. The data is consistently above 
that predicted by McAdams equation and can be represented by: 

Nu » 0.017 Re 0,8 ^ Pr 0 * 1 * (2) 

Subcooled nucleate boiling data is plotted in Figure 6. Figure 6 

contains data at the same flow rates, P — 60 psia, L/D ■ 37 ,h> and a 

range of subcooling from about 0 to 175°F. Three other curves are 

plotted to aid in describing the results. 

First, the curve of (q/A)pQ is the curve of heat flux versus wall 

superheat (T -T ) obtained from non-boiling prediction. This curve 
w sat 

is based on G = 1.19 x 10^ lb/hr ft 2 , and T su ^ « 109-183°F, where 



Numbers denote reference listed in Bibliography 
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( pc “ h <VV “ h ^ T w" T sat + T sat “ T b^ and h is determined Iron 

properties at bulk temperature corresponding to the data using Eq (2), 
Second, the curve of (q/A)^ represents the locus of points at which 
surface boiling begins and where the heat flux departs from the non- 
boiling prediction. This is discussed in Reference (5) and is given by: 



(q/A) i 



15.60 



1.156 

P 




T sat> 



2.30/p 



0.023U 



(3) 



The curve plotted here is for p » 60 psia. 

The third curve represents the mean line through data of Vessel 
for fully-developed subcooled nucleate boiling and vias determined as 
indicated in Figure (7). 

The data covers a range of flow conditions in the transition region 
from the inception of boiling to a point on the fully-developed boil- 
ing curve at which bulk boiling begins. As fully-developed boiling is 
reached, the data approaches an asymptote which agrees with Vessel's 
data. Figure (6) illustrates that in the transition region heat trans- 
fer rates are dependent on both flow rate and subcooling, and as the 
fully-developed boiling region is approached the dependence on flow 
rate and subcooling disappears. 

Incipient boiling does appear to begin at the point predicted by 
the intersection of the incipient boiling curve and the non-boiling 
prediction. 

Bulk boiling data is presented in Figures 8 and 9* Figure 8 shows 
data for L/D - 37. U, at the two lower flow rates for bulk boiling 



(i.e., at points where the bulk fluid was nt the saturation tempera- 
ture). The data represents the range of steam quality from 0 to 2k%. 
The fully-developed curve is again plotted for comparison. The data 
of Figure 8 to appear to follow the fully-developed boiling curve for 
the range of quality indicated. 

Figure 9 contains data obtained with T.S. -2, 3, and 3 modified. 

The range of steam quality is extended to almost U0$, and data is 
plotted for varying L/D r s. The data for a particular L/D can still 
be represented by a form of Eq. (1). But, as quality increases (rep- 
resented by increasing L/D), the fully-developed boiling curves are 
shifted toward the lower wall superheat. The effect of quality seems 
to override the pressure effect, i.o., the lower pressure results in 
higher wall superheat. 

Similar results appear in Figure 10 for data with subcooled 
boiling conditions. As L/D increases, at first, wall superheat in- 
creases, and then decreases. This appears to be independent of the 
degree of subcooling within the range covered in Figure 10 ( T g * 

30 - l50°F) o 

In an attempt to determine the effect of steam quality on heat 
transfer coefficient, Figures 11 through lU show the experimentally de- 
termined heat transfer coefficient versus steam quality. Figure 11 con- 
tains data obtained with T.S. -2 for G ■ 0.59 and 1.19 x 10^ lb^/hr ft^. 
The results indicate that heat transfer coefficients are dependent on 
heat flux, and flow rate alone in the range of steam quality. At 
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quantities above 10-15$ the data appears to merge for each flow rate 
and the dependence upon heat flux disappears. The effect of flow rate 
is quite pronounced as would be expected. At lower flow rates the effect 
of quality on the local transfer coefficient becomes apparent at lower 
values of quality. 

Figure 11 shows the results of T.S. -2 and 3 at G « 1.19 x 10 6 , and 
Figure 12 shows similar results for G - 0.59 x 10 6 . Data obtained from 
these test sections appears to be inconsistent. Results of T.S. -2 in- 
dicate an increase in h as flow proceeds downstream. If the data ob- 
tained from T.S. -3 is then compared with this, it would appear that h 
is relatively constant between 37 . h L/D and 57.9 L/D. Beyond that point 
h reaches a maximum and then decreases. Because of the inconsistency in 
the data, T.S. -3 was modified. The thermocouples located upstream on 
T.S. -3 (L/D - 57.9 and 6U.1) were moved to positions which corresponded 
to the downstream thermocouple positions of T.S. -2 (L/D - 30.9 and 37. h). 

Data taken with this test section is shown in Figure 13. Also plot- 
ted for comparision is data taken with T.S. -2 and 3 at G - 1.19 x 10^ 
lbs/hr ft . Data for each test section was taken at approximately equal 
values of heat flux and inlet temperature, in so far as data was avail- 
able from the other test sections. Results of these last tests compare 
favorably with that obtained from test sections 2 and 3. The absolute 
values of h at corresponding L/D's are about the same. The only differ- 
ence is that data from the modified test section shows an increase in h 
at the last thermocouple position instead of the decrease noted with 
T.S. -3 
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In order to compare the bulk boiling data with previous work in 
this area, Figure lli shows the dimensionless ratio h/h^ versus the 
Martinelli parameter (l/X^^) , where the Martinelli parameter is defined 
by: 

• 1/X U - (x/l-x)°-9 (y> L /y? v )°* 5 (U) 

h is the experimental heat transfer coefficient and l^is the liquid heat 
transfer coefficient as given by the McAdams correlation. Also plotted 
here is Chen's F function and his prediction of h at parameters approx- 
imating the mean values for the data. The curve lies below the data at 
the lower L/D's of T.S. -2 , and 3 modified, but above the data for high 
L/D's of T.S. -3 and 3 modified. Chen's correlation is presented in a 
later section. 

Figure 15 shows similar data at G = .59 x 10^ along with Chen's F 
function and his predition of h at parameters approximating the mean 
values for the data. Here Chen's F factor seems to divide the data for 
low L/D's and high L/D's. At the low flow rate Chen's prediction of h 
is higher than the experimental values. 
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DISCUSSION OF RESULTS 

Critical Heat Flux 

In the course of obtaining this data several interesting phenomena 
occurred. Although the primary intent of this study was to obtain heat 
transfer data, it is considered necessary to mention these items briefly. 

First, the fact that (qA) cr increased as much as 50$ as the exit 

quality increased from 0 to about 25$. Previous studies in this area 

with subcooled fluid conditions at the exit have indicated that (q/A) 

decreases as the amount of exit subcooling decreases. The most recent 

studies show that as zero subcooling is approached (q/A) appears to 

cr 

level off and in some instances increases . This is indicated by the 
data of VJessel shown in Figure U. One of the best known studies of 
burnout with exit quality is that of Lowdermilk (6). His results in- 
dicate that (q/A) cr is at a maximum when exit quality is approximately 
20$ depending on other flow parameters as well. The data did not cover 
the low range of exit quality. Hence, previous work indicated that, at 
a constant pressure and flow rate, ( q/A) decreased as the exit subcool- 

ing decreased. Critical heat flux reaches a relative minimum, increases 
in the low range of exit quality to a relative maximum, and then falls 
off sharply. Bergles (l) also shows an increase in critical heat flux 
with an increase in quality up to 12$. Hence, the data presented here 
is in agreement with other studies. This critical heat flux is not as 
high as that attained with subcooled exit conditions. For more infor- 
mation in this area the reader is referred to recent studies by Loosmore 
and Skinner (7) on critical heat flux with subcooled boiling and by 
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Lopina (&) on critical heat flux with exit quality conditions. 

Second, the actual location of burnout of the test section was 
not always at the exit. In fact, at exit quality near zero, and with 
exit quality about 25 %, burnout did occur at the exit. Between these 
limits, the location of burnout was considerably upstream of the exit. 

In one run with exit quality apnroximately 13.6a burnout occurred at 
18.7 L/D, while the overall tube dimensions in all runs was Iil.U L/D* 
Calculations indicate that the bulk fluid condition at the burnout 
location was slightly subcooled. This indicates that in the low range 
of exit qualities, exit conditions may not be the correct parameter to 
consider in determining critical heat flux. Lopina (8) and Waters, etal 
(9) have also observed upstream burnout and discuss this in more detail. 
Waters* study shows that upstream burnout is reproduceable and not neces 
sarily caused by instability. 

Third, it was observed during operation that as the exit quality 
increased, the pressure drop from the test section exit to the exit 
plenum chamber became quite large. Hence, when measuring critical 
heat flux with net vapor generation, exit pressure must be determined 
at the test section exit and not in an exit plenum chamber. Results 
of this study indicated that there was essentially no pressure dif- 
ference between these two locations at zero quality, but at 17.8a 
quality the pressure difference was approximately 12.5 psia at G - 
1.19 x 10^ lbs/hr ft^. To explain this it is necessary to give a 
brief discription of the method used to control the pressure. 

The exit pressure level was controlled by means of a plug valve 
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downstream of the exit. At low values of heat flux, this valve was 
about midway between the fully opened and closed positions in order 
to maintain the desired pressure. As the heat flux was increased, 
and boiling conditions were reached, the exit pressure began to in- 
crease. This valve was then opened to reduce the pressure to the 
desired level. When bulk boiling conditions were reached, it was 
observed that pressure control was no longer possible, since the 
valve was wide open. As indicated previously, when the longer test 
section was used, pressure control was lost prior to obtaining the 
maximum heat flux for that series of runs. This then illustrates the 
choking effect that causes the pressure level to increase in the test 
section. 

As the steam-water mixture leaves the test section, it enters a 
plenum chamber made up of 3/U inch brass pipe. Expansion of the steam 
can take place immediately. The pressure drop encountered from the 
test section exit and the plenum chamber is the result of this rapid 
expansion. 

Heat Transfer Results 

Non-boiling data shown in Figure 5 plots consistently above the 
correlation of McAdams. The divergence appears to be primarily due to 
the radial property variation which is not properly accounted for in 
the correlation. This temperature difference effect is clearly demon- 
strated in Reference (10). In any case, Eq. (1) is not intended to be 
a generally valid correlation, but is merely indicated to represent the 
data. 

Figure 6 shows the portion of the boiling curve from the inception 
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of boiling to the fully-developed boiling curve for three flow rates 
and two ranges of subcoolings. As indicated previously, inception of 
boiling did occur at the predicted value of wall superheat and heat 
flux. Beyond this DOint, the slope of the curve increases until fully— 
developed-boiling is reached. Using the data of Vessel to establish 
the fully-developed boiling curve (Figure 7), the data does approach 
this line as an asymptote. Kence, the subcooled data appears to be con- 
sistent with the work of previous studies. 

At the onset of this study. Professor Bergles pointed out, that at 
high heat fluxes with low subcooling, the data appeared to be shifted 
from the fully-developed curve, toward lower wall superheats. It has 
already been mentioned, that fully-developed subcooled boiling data 
(Figure 10) show a change in wall superheat at different positions a- 
long the length of the tube. At the inlet, subcooled liquid velocity 
can be specified as V]_. As the fluid moves down the heated sections, 
surface boiling begins and non-equilibrium voids are formed. Therefore, 
the volume occuoied by the liquid is reduced. In order that continuity 
be maintained, the liquid velocity must increase. Let this new velocity 
be \ r 2 » As boiling becomes more vigorous along the tube, more volume is 
occupied by vapor and the liquid velocity continues to increase. This 
illustrates how the liquid velocity may increase along the length of the 
tube. With this in mind, and noting from Figure 6 that by increasing 
fluid velocity, fully-developed boiling is reached at higher values of 
heat flux. 

If we consider three fluid velocities which do not differ by much, 
and determine the heat flux at which the lowest velocity reaches the 
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f ully-de velo ped boiling curve, it is clear that at this same heat 
flux, the wall superheat will be reduced at the two higher velocities. 

This example shows how the wall superheat can decrease at higher L/D's 
with subcooled liquids. This also indicates that the local heat trans- 
fer coefficient must increase in this range. The heat transfer co- 
efficient is defined as: 

h - (q/A) / (T w - T b ) (5) 

As long as the bulk fluid is sub^ooled, the bulk temperature is in- 
creasing along the direction of flow. If (T w - T sat ) decreases as 
indicated previously and T sat decreases with pressure, then T w must 
also be decreasing. As a result, (T w - T fe ) must decrease and the heat 
transfer coefficient increases due to the presence of void fractions. 

This explains the trend toward lower wall superheat in data at high 
L/D's as shown in Figure 10. 

Bulk boiling data (Figure 9) also shovjs a reduction in wall super- 
heat with increasing L/D. Bennett, etal. (11) presents a discussion on 
the effect of net vapor generation, similar to the discussion above on 
the effects of void fraction on fluid velocity. This explains the lower 
wall superheats obtained at the exit of the test section. In the region 
of bulk boiling where the bulk temperature is the saturation temperature, 
a decrease in wall superheat means an increase in h for a specified heat 
flux. Therefore, wherever net vapor exists, there is an increase in the 
local heat transfer coefficient. 

Now, consider a section along the tube where the bulk fluid is slight- 
ly subcooled. As flow proceeds down the tube the bulk fluid approaches 
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saturation conditions. The rate at which bubbles nucleate at the 
surface decreases since the fluid velocity is greater. At the same 
time, the bubbles, which have just departed the surface, do not 
collapse and the volume fraction of vapor remains nearly constant. 

Then, the fluid velocity will also remain constant and the degree 
of wall superheat predicted by the boiling curve will not change. 
Therefore, in this region near zero quality, the heat transfer co- 
efficient is relatively constant as well. Further along the heated 
section, net vapor generation continues by means of surface boiling 
and/or evaporation at the liquid vapor interface. At this point, the 
fluid velocity begins to increase again and continues to increase 
until the test section exit is reached. This appears to explain the 
data of Figure 9> as well as that of Figure 11. 

Figure 12 with data of T.S. -2 and 3 at G - 1.19 x 10 6 lbs/hr ft 2 
shows the discrepancy in the data obtained from these two test sections. 
Data obtained with T.S. -3 shows that, as x increases, h reaches a max- 
imum and then decreases. This could not be explained adequately. Con- 
sidering data for an approximate equal heat flux for the two test 
sections, there also appears to be a discontinuity in the data between 
L/B * 39.U and 57.9. The value of h increases as L/D increases to 39.h, 
but apparently there is no further increase in h to L/D m 57.9. Beyond 
this point, h reaches its maximum value and decreases. Results for G » 
5.9 x 10^ lbs /hr ft 2 shown in Figure 12 are similar. 

It was expected that data obtained with T.S. -3 modified would re- 
solve this discrepancy, figure 13 shows data for one flow rate and one 
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value or heat flux for T.S. -2, 3> and 3 modified. The results of this 
show that the absolute magnitudes of h remained about the same. The only 
difference which did appear consistently was that at high L/D's, h in- 
creased rather than decreased. The heat transfer coefficient then show- 
ed and increase over the first half of the T.S. , then remained a rela- 
tively constant value, and finally increased again near the exit. One 
possible explanation is; at the lower L/D's flow conditions had juSt 
reached the point at which the void fraction is nearly constant. Over 
a short range beyond this point, h is constant with x. If this were 
true, the results as indicated in Figure 11, are misleading and quality 
does not significantly affect h until it reaches a value of about 20%. 

Another method of evaluating the data would be to draw the mean 
line through the data from both test section. It is uncertain which 
data, if any, is actually incorrect. If this were done in Figure 12, 
it would indicate the effect of quality becomes apparent at about 15$, 
but its overall effect would not be as great as indicated in Figure 11. 

Figures 1$ and 16 show the same data on another set of coordinates, 
h/h^ versus l/X^,„. Dengler and Addoms (12) were the first to use these 
coordinates in evaluating this type of data. Here h^ is defined as the 
conventional liquid heat transfer coefficient determined by McAdams* 
correlation. They proposed a correlation for two phase heat transfer 
coefficients in terms of l/X (the Martinelli parameter) and an F 

1 x 

factor to account for the degree of nucleate boiling present. As the 
fluid progresses into the annular flow regime, with continuiously in- 
creasing steam quality, the increased velocity of the two phase mixture 
induced by the vaporization process suppresses the nucleate boiling proc- 
ess. Beyond this point the heat transfer coefficient becomes governed 
by the forced convection process. 
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liitih "this in mind Chen (13) proposed the following correlation 
for the two phase heat transfer coefficient: 



^ “ Wiic + ^mac 



( 6 ) 



where h mic is the microconvective portion due to the boiling mech- 
anism and hence dependent on liquid and bubble parameters: 



\ic ■ - 00122 



^ g 0 °-« 



r°'Ul' 29 \°- 2h /> v °- 2k 



x (AT sat )°- 2U (AP sat ) 0 - 75 s 



(7) 



and h is the macroconvective given by: 

III ext* * 



h mac - h, F (») 

Again h^ is given by McAdams’ equation. F is a function of the 
Kartinelli parameter and is defined as the ratio of the two phase 
Reynolds number to the liquid Reynolds number. S is called the 
suppression factor and approaches unity at zero flow rate, and zero 
at infinite flow rate. Gouse (lii) suggests that this suppression 
factor should go to zero much sooner and thus acknowledge the possi- 
bility of total suppression of nucleation at finite Reynolds number. 
Both S and F are presented graphically in Reference (13) » 

Chen restricts the applicability of this correlation to systems 



with annular or annular mist flow in vertical geometries. This 
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correlation successfully predicts tho data of numerous independent 
studies to within plus or minus l5/£, and seems to be one of the best 
available. Gouse (lit) also states that there is no reason why it 
should not apply to annular flow in horizontal tubes. 

Data obtained with thermocouples at the top and bottom of the 
tube wall, at the same location along the length, produced no indica- 
tion of stratification. Thermocouples indicated that the temperature 
difference obtained in all cases was less than the accuracy of the 
measurement. Also, the greater temperature fluctuated between the top 
and bottom. With this in mind Chen's correlation is considered a 
suitable means of checking the data. 

As noted earlier Chen's prediction seems to plot between the data 
for high and lowL/D's indicating that the difference can be considered 
as scatter in the data. It must be noted, that scatter of such a con- 
sistent nature is unusualo 
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CONCLUSIONS 

Critical Heat Flux: 

The critical heat flux increased as much as $ 0 % as exit quality 
was increased from 0 to 25$. 

Heat Transfer 

1. The local heat transfer coefficient, obtained at subcooled 
boiling conditions, does appear to increase due to the pre- 
sence of voids formed during the boiling process. 

2. As bulk fluid approaches saturation conditions the local heat 
transfer coefficient is constant over a low range of steam 
quality. The extent of this range increases as the flow rate 
increases. 

3. For the flow rates considered in this investigation, the local 
heat transfer coefficient shows a significant increase as qual- 
ity is increased beyond 15$. 

U. Chen’s correlation appears to be iairiy accurate in the low 
range of steam quality, but it's accuracy improves as quality 
increases and nucleate boiling is suppressed. 



-27- 



RECOMMENDATIONS 

Critical Heat Flux 

1. "When bulk fluid conditions at the exit contain vapor, the 
exit conditions must be determined by the pressure at the 
test section exit rather than the pressure in an exit plenum 
chamber. 

2. Since upstream burnout is possible with exit quality conditions, 
and stable flow, it may be of interest to determine if some para- 
meter, other than exit quality might give a more meaningful in- 
dication of the results. 

Heat Transfer 

1. Additional experiments should be conducted with test sections 
similar to T.S. -3 and instrumentation covering a greater por- 
tion of the test section. This is suggested to see if there 
is actually a region along the test section. where the local 
heat transfer coefficient is constant. 

2. Additional study of heat transfer coefficient in the low qual- 
ity range is still required, in order to develope a more accu- 
rate prediction in the range where heat transfer rates are de- 
pendent on the combined effects of boiling and convective heat 



transfer. 
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APPENDIX A 

EXPERIMENTAL APPARATUS 
Flow Loop and Power Supply 

The flow loop used in the experimental program was one already in 
operation at the M.I.T. Keat Transfer Laboratory. This facility is 
shown schematically in Figure 1. It was a closed loop system in which 
all components are of corrosion resistant materials. The system con- 
tained the main circulating pump, an accumulator, the test section line 
with parallel flow meters, preheater, test section and its instrumenta- 
tion, the bypass line, and a heat exchanger utilizing city water. Aux- 
iliary equipment included a fill pump, supply tank, degassing tank, and 
a continuous demineralizer. The test section power was supplied by an 
a.c. motor - d.c. generator set. A discription of the major components 
follows. More detailed information is given in Reference (1)« 

The main circulating pump was a two-stage, turbine-type pump driven 
by a 3hp induction motor. It provided a head of 250 psi at 3.6 gpm. A 
bladder type accumulator, pressurized with nitrogen, located at the pump 
exit served to damp out pressure fluctuations. The ball valve in the by- 
pass line was used to control the by-pass flow rate and pressure, which 
in turn determined the pressure in the test-section line. The test sec- 
tion line contained a Fischer-Porter flow meter and a set of preheaters 
in series. There are four 5kw preheaters plus an additional 5kw emersion 
unit fabricated in a 3 inch copoer tube and provided with 0 to 100$ con- 
trol by a Powers tat auto-transformer. This arrangement allowed varying 
preheat from 0 to 25kw, one being variable and the other allowing incre- 
ments of 5kw each. 
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ihe test, section power was supplied by motor— generators. The 
generators were driven by hhO volt, 3 phase synchronous motors. Two 
36kw d.c. generators, each nominally rated at 12 volts and 3000 amps 
were connected in series. The generators -.'ere provided with water 
cooled shunts in parallel with the test section which allowed them to 
be run open circuited at the test section for starting, or after burn- 
out. 

Instrumentation 

Instrumentation was available for reading pressure levels, differ- 
ential pressure, temperature, test section flow rate and test section 
voltage and current. Pressure levels were read on bourdon- tube gages 
as shown in Figure 1. At the test section inlet and exit were 200 and 
100 psi test gages ’with an accuracy of plus or minus l/k% of full scale. 

Test section flow rate was determined from the larger flowrater meter 
which had interchange ole tubes and floats. The ranges covered are 20- 
160, 35-UOO, and lU0-13o0 lb/hr. All units were calibrated as installed 
in the system. 

The power supplied to the test section was obtained from the test 
section voltage and current. The voltage was read directly on a Weston 
multiple-range d.c. voltmeter with a specified accuracy of plus or minus 
l/2%. The current was determined from a standard shunt with a calibra- 
tion of 60.17 amp/mv. 

All temperatures ’were measured with copper-constantan thermocouples 
made from 30 gage duplex wire. The fluid bulk temperature at the inlet 



•and exit of the test section was measured by inserting thermocouples 
directly into the fluid using conax fittings with lava sealants. 

Methods of measuring wall temperature are deecri'oed later. The out- 
put. voltages of the thermocouples and the shunt were displayed on a 
continuous recorder. The recorder was a Brown, single channel in- 
strument having ranges of 0-6, 5-11, 10-16, 15-21, and 20-26 mv. 

The manometer system used to measure differential pressure con- 
sists of two Meriam 60 inch U-tube manometers, manifolds and valves, 
and connecting lines of rubber hose. The system is designed to read 
a maximum of 10 pressure differences with either one or two reference 
points on either manometer. One manometer was filled with mercury 
which gave a maximum range of ‘approximately 25 psi. The other mano- 
meter contained an oil with a specific gravity of 2.00 which allowed 
a maximum cressure difference of approximately 2 psi. A typical test 
section with three pressure taps is shown in Figure 2. The pressure 
drop from the reference "A" located at the inlet, to points 2, and 3 
could be read on either the mercury or oil manometers depending on the 
magnitude of the pressure drop. The two pressure gages were connected 
to the inlet and exit to give the pressure level and a check on the 
pressure drop. Vents were included so that all lines could be purged 
of air. 

Measurement of Tube—Vali Temperature 

The most difficult problem in any heat transfer experiment is the 
measurement of the temperature of the heat transfer surface. The basic 
method of measuring tube wall temperature is described in Reference (1)„ 
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The method used here was as follows: 

Heated shields constructed of standard aluminum tubing, with 
threadlike grooves cut into the outer wall, were wrapped with in- 
sulated Chromel A-26 wire. Spacers fabricated from aluminum served 
to hold the shield in place. 

Power was suppled in the Chromel heating coil by a 110 VAC supply 
which is stepped doxm to about 25 volts by a Stancor filament trans- 
former. The power to each shield is controlled by a Powerstat. This 
arrangement permits fine control of shield heating. 

A layer of Scotch Electrical tape Number 2? was wrapped around 
the tube at each thermocouple location. The measuring junction was 
then bound radially around the tube with another layer of tape. Thermo- 
couples were also taped to the inside of the shield. Sufficient length 
of all thermocouples were left inside the shield to avoid conduction 
errors. The space between the shield and tube was filled with asbestos 
fiber. Precautions were taken to insure the shield was electrically 
insulated from the test section. 

Power to the shields is adjusted until the shield thermocouples 
and the tube wall thermocouples within the shield section indicate 
approximately equal Eli?. 

Test Section 

The test section was constructed from 30U stainless steel tubing 
with an inside diameter of 0.2U25 inch. Brass bushings were cut from 
3/U inch round stock. Holes were drilled for the test section. The 
bushings were then sweated to the tube with silver solder. Holes 0.099 
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inch were then drilled in the bushings for pressure taps. These 
holes did not penetrate the test section. A number 75 drill (0.021 
inch) was used to drill the pressure tap hole. Beryllium copper pres- 
sure tap tubes (2mm ID) were then set in the bushings and soft soldered 
m place. The inside of this tube was filled with milk of magnesia to 
prevent the tap being plugged. With this completed, the milk of magnesia 
could then be removed. 

An additional pressure tap was placed near the center of the test 
section. Here a short stainless steel tube was first filled with milk 
of magnesia and then silver brazed to the test section. Excess braze 
was then filed away to give a very small fillet at the test section. 

The milk of magnesia was removed and the pressure tap hole drilled (0.021 
inch) through the test section. Any burrs inside the tube were removed 
with fine emery cloth and pieces of steel wool pushed through the tube. 

Thermocouples were then attached to the test-section wall as des- 
cribed previously. Thermocouples were located at four positions along 
the length. At two locations the thermocouples were placed perpendicu- 
lar to the tube axis, and the lead wires wrapped around the tube and 

taped in place. At the other two locations, thermocouples were again 
placed perpendicular to the tube axis, but at each of these positions 
two thermocouples were used, one on the top of the tube and the other 
at the bottom. This set up was intended to detect stratification. 

The guard heaters were then assembled onto the test section. After 
the first guard heater was in place a section of 0.099 inch OD tubing 

was placed over the center pressure tap and soft soldered to the stain- 

less steel tube to extend this pressure tap. A diagram of this test 
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section is shown in Figure 2. 

T.S. -1 and 2 were identical in construction. Accidental burnout 
early in this investigation is the reason for this. The location of 
pressure taps and wall thermocouples is indicated below. 



T.S. 1 and 2 
Gaining length 6 in. 

Heated length 10 in. 

Center pressure tap at 27.3 L/B 



T.S. 3 

Calming length 6 in. 

Heated length 20 in. 

Center pressure tap at li5.5 L/D 



Distance to thermocouples 
No. 1 lli.U L/D 

2 23.7 

3 30.9 

h 39. U 



1 57.9 L/D 

2 6U.1 

3 72.ii 

ii 78.5 



T.S. 3 modified 

Calming length 6 in. L 

Heated length 20 in. 2 

Center pressure tap at li5 • 5 L/D 3 

h 



30.9 L/D 
39.U 
72.ii 
78.5 
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APP5K3IX I 
SAMPLE CALCULATION 



1 . Calculation oi' heat input to the test section; 

(a) from the current and voltage measurements 

q„ “ 3.U13 SI - 3.H3 (60.17) E E- 

= 203.36 E E g 3tu/hr 

S and E were measured value s 0 
s 



(b) from current and resistance 

q K 3 3«ul3 I 2 R 1.236 x 10^ S 2 R 

e s 

R = ^ L/A c 

For T.S. -1 and 2 L = 10/12 ft, 

A c = (If/ 4 ) (ID 2 - 0D 2 )/lUu - 2.113 x 10^ ft 2 

R “ 0.39U X 10^ 



liras determined from Reference (1) as a function of the 
tube temperature. The tube temperature is given by; 



T t ■ T o W - <^ T »/2) 

\T was determined from Figure 2 and is a function of T 
w ow 

I. T qw was measured by means of thermocouples converted with 
N.B.S. tables. 



(c) from mass flow rate and bulk temperature change 
<3 C w ^p ^exit - "i^ 

T ex it and as well as all other temperatures was obtained 

as indicated above for T . 

o w 

A sample set of data and calculations are shown on the following pages. 
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3. Average heat flux q/A = I;80 93/ .0529 M 9.159 x 10^ Btu/hr ft 2 
U. The change in enthalpy through the test section is 



Ah’ - q h - 135 Btu/lb 



m 



5. The change in bulk temperature from inlet to any position along the 



T,' i3 T^ 3 Tq and Tyj denote guard heater temperatures. 

The a and b indicate temperatures at the top and bottom of the tube 
respectively. 



fl 

















tube is given by: 
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A b “ (^exit “ ^i) x '/ ^ 

6. The pressure at any point along the tube was obtained graphically 
as indicated earlier. 

7. Entering the steam tables with the local pressure, the saturation 
temperature, the saturated liquid enthalpy ar.d the latent heat of 
vaporization were obtained. 

8. All fluid properties at each thermocouple position were obtained 
at the local bulk temperature. These include Cp, k and P r# 

9. At each position the following values were obtained as indicated 



above 


• 










L/D 


x»A 




m ** 

x ow 


AT 

V 


<*n 
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.35 


a 
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23.7 
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565.7 






37. 5 


.906 


b 


563.0 


115 


359.3 



P T 
x ’ sat , 

x' 

60.2 292.9 

59.9 292.6 

59.6 292.3 

5<?.5 292.1 



The temperatures at a and b were averaged to give the outside wall 
temperature at these locations. 
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L/D 


a m 
^ x b 
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10. The local value of h was not calc 
but was determined for non-boilin- 

h - (q/A) / (T w - 

11. Nu and Re were determined for non- 
Nu/?r C °^ . 



A T . 
suo 


< T w - T b ) 


191.6 


21:2.7 


160.6 


215.0 


U:7.3 


200.0 


126.2 


163. h 



elated for subcoolcd boiling data, 
; and bulk boiling data as shown: 

V 

-boiling data as well as the ratio 



Nu = hO/k - (h/k) (0.0202) 

Re - GD/m = 2.Ji0ii x 10^ /jU[ ©• w - 361: lb m / hr 

Pr = C p ^/k 

Pr is tabulated as a function of temperature in Reference (15) with 
the other properties for water. 

12. Bulk boiling data included some additional calculations as indicated: 

(a) local enthalpy h r x t = h*^ +^h' v| 

where Ah' . = (h* - h l . ) x’/l 

NOTE: When this value exceeded the value of saturated liquid 

enthalpy, net vapor was present. The local quality was 
calculated using the local enthalpy, saturated liquid 
enthalpy and latent heat of vaporization. 

(b) Evaluation of the Kartinelli parameter required values of JA, 
and p for steam. Mv was obtained from Reference (l6) and 
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was obtained from the steam tables. The Martinelli 
parameter was then calculated for each data point at which 
net vapor was present. The Martinelli parameter is defined 
as: 

v% - (-2-) 0 * 9 (J= L /P V )° A 

(c) To obtain the ratio h/h^, the liquid heat transfer coefficient 
(h^) was calculated at each of these data points with McAdams 
equation: 

h L - 0.023 (k/D) Re 0 " 8 Pr 0 *^ 

This equation was used because it provided a means of compar- 
ing the data with Chen's '.cork in this field. His value of h^- 
is equal to h^, as determined above, multiplied by his F factor. 



A summary of the data is contained in Appendix C 
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APPBNDIX C 
STJKKARY OF DATA 
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FIG. 8 BULK BOILING DATA T.S. 2 
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